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q ■ Abstract 

In determining the mesonic wave function from QCD inspired potential model, if 
■ the linear confinement term is taken as parent ( with columbic term as perturbation ), 

P^H ■ Airy's function appears in the resultant wave function - which is an infinite series. In 

| the study of Isgur-Wise function (IWF) and its derivatives with such a wave function, the 

infinite upper limit of integration gives rise to divergence. In this paper, we have proposed 
some reasonable cut-off values for the upper limit of such integrations and studied the 
subsequent effect on the results. We also study the sensitivity of the order of polynomial 
approximation of the infinite Airy series in calculating the derivatives of IWF. 
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^h! 1 Introduction: 

O; 

The quark composite systems with one heavy quark have been the focus of interest for the 
last few decades. In infinite heavy quark mass limit, the heavy quark sector of QCD becomes 
CNl independent of quark masses and the effective Lagrangian of the heavy quark effective theory 
( HQET ) exhibits additional spin flavour symmetries and this simplifies the calculations of 
matrix elements of electroweak transitions. The semi-leptonic transitions between heavy mesons 
^ ■ are of great importance. In the limit of infinite quark masses, all the mesonic form factors can be 
expressed in terms of a single universal function, called the Isgur-Wise Function [1]. This IWF, 
the universal function describing a class of semi-leptonic transitions, has been studied for many 
years. The shape (or form) of the IWF and its derivatives (slope and curvature) at zero recoil, is 
essential for determination of CKM matrix elements and that requires a reasonable description 
of the IWF. The main part of the IW function is the wave function of the hadrons. The wave 
function for the heavy-light mesons have been calculated earlier within the framework of QCD 
potential model (Cornell potential ) [2,3] with considerable accuracy [4]. This has been deduced 
both with columbic term in potential as parent [5] and also with the linear confinement term 
as parent [6]. The characteristics like slope (charge radii ) and curvature (convexity parameter 
) of IWF have been reported in both the two cases with certain limitations. In the present 
paper, we have considered the linear confinement term (br) as parent and Columbic term as 
perturbation. This has been done applying Dalgarno's method of perturbation theory up to 
first order correction [7]. As Airy's function appears in the wave function for unperturbed 
Hamiltonian, the integrability of the otherwise divergent infinite Airy's function series in the 
IWF is a question of consideration in the present work. We have introduced some reasonable 
cut-off for upper limit of integration of IWF keeping in mind the nature of Airy's function 
and the boundary condition of IWF ( £(!)= 1) and studied the variation of results with this 
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cut-off value at different orders of Airy's function. We also study the sensitivity of the order 
of polynomial approximation of the infinite Airy's function when compared with experimental 
result of the derivatives of IWF. The section 2 contains the essential formalism, section 3 , the 
calculations and results and the section 4 contains conclusion. 



2 Formalism: 

2.1 Potential Model: 

The potential under consideration ( Cornell potential ) is : 



V (r) = -^ + br + c (1) 
3r 

Here we take br as parent so that our unperturbed Hamiltonian [8] is 

* = + * (2) 

with 

ff'^ + c (3) 
as perturbation. Here \i is the reduced mass, which is 

H= (4) 
m q + mQ 

We take the value of b to be 0.183 GeV 2 from charomonium spectroscopy [9] and constant C 
to be 1 GeV . It is to be mentioned that, In the infinite heavy quark mass limit ( mg — > oo ) , 

a = lim — — m„ (5) 

mg^oo m<? _|_ mQ 

Under this consideration, the two body Schrodinger equation [11, 12] for the Hamiltonian 
H = H + H' is : 

H\V >= (H + H')\V >= E\V > (6) 

2.2 Wave Function: 

To find the unperturbed wave function corresponding to Hq we get the radial equation for 
potential br for ground state S (1 = 0), following the formalism of reference [13] , as : 

Id 2 2d 

[- — (-T^ + --T-) + ^r] = ER(r) (7) 
L 2fi K dr 2 rdr J 1 y ' y ' 

where R(r) is the radial wave function. We introduce u(r) = r R ( r ) and dimensionless variable 
p , where - 

p = (2pb)^r - (fy/*E (8) 
The above equation then reduces to: 



The solution of this second order homogeneous differential equation [f 4] contains linear combi- 
nation of two types of Airy's functions Ai[r] and Bi[r]. The nature of the Airy's function [15,16] 
reveals that 



Ai[r] — > as r — > oo and 
Bi[r] — )> oo as r — > oo. 

So, it is reasonable to reject the Bi[r] part of the solution. The radial wave function has thus 
the form: 

u{r) = NA l [{2pb) 1 '\r-^)} (10) 

where N is our normalization constant which has the dimension of GeV 1 ^ 2 .The boundary 
condition u(0) = [17] gives us the unperturbed energy for ground state [13]: 

W^ = E=-(^-f 3 Po (11) 

where p is the zero of the Ai[r] , such that Ai[p ]=0 [18]. 
p has the explicit form - 

37r(4n — l) l2 /-? . „ , 

po = —[ J ' , m our casen = 1 j or ground state. (12) 

From this we get the unperturbed wave function for ground state to be : 

^{r) = -^A l [(2ub)^r + p } (13) 

^(r) = ^=-A l [p 1 r + Po ] (14) 
where we have taken Pl = i^lub) 1 ^. 

The first order perturbed eigen function vf^ 1 ) and eigen energy W 1 can be calculated using the 
following relation: 



Ho^' + H'V = W {0) y' + W {1) ^° (15) 



We find, 



W m = / r 2 H' | ^° | 2 dr (16) 
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Employing Dalgarno's method [19] , the first order wave function comes out to be [6]: 

^(i) (r) = _!^^o + 6r + c) (17) 

where a , a>\ and a 2 are terms which involve a s , b, p, E and c. These are having dimensions 
of GeV 1 / 2 , GeV 3 / 2 , GeV 5 ^ 2 respectively and have explicit form, considering Airy order up to r 3 , 
as given below [6]. 
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0.8808 (V) 1/3 «2 4W /1 x 0.21005 

a ° ~ (E-c) fi(E - c) + 3a s (E-c) (18) 
ba 4W 1 x 0.8808(bfiy/ 3 0.6535 x (bp) 2 / 3 
ai ~ {E-c) + 3a s {E - c) (iT^) ( ™ 

x 0.1183 

a 2 = — ~ (20) 

Sets 

The total wave function with first order correction is: 

*tot(r) = *°(r) + V 1 (r) (21) 

which upon substitution yields- 

%ot(r) = N 'i^M(Pir + Po ]-^(^ + br + c)] (22) 

Here N' is the normalization constant of total wave function which is also having the dimension 
of GeV 1 / 2 . Considering relativistic effect on the wave function, the total relativistic wave 
function is given by [11]: 

*rej(r) = * tot (r)(-)- £ (23) 



This gives, 
Here, 



<Mr) = N'i-^A^r + p ] - ^(^ + br + c)](-)" e (24) 



a »=zk ande=1 -f- ( f )2 (25) 



2.3 Isgur-Wise Function: 

In case of semi-leptonic decay of hadrons ( mesons ) , in the infinite mass limit, a new symmetry 
called spin-flavored symmetry, will emerge and the Heavy Quark Effective Theory ( HQET ) 
will be suitable. In this theory, the strong interactions of the heavy quarks are independent 
of its spin and mass [20] and all the form factors are completely determined, at all momentum 
transfers, in terms of only one elastic form factor function, the universal Isgur-Wise function 
£(v,v'). £(v,v') depends only upon the four velocities v v and v v i of heavy particle before and 
after decay. This £(v,i/) is normalized at zero recoil [21]. If y = v v .v v > , then, for zero recoil 
(y=l), £(y) = 1. With increasing recoil y grows. In explicit form IW function can be expressed 
as : 

i{y) = i- P 2 {y-i) + C{y-if + (26) 

p 2 is the slope parameter at y=l given by - 

P=— Ejr \y=l (27) 

p is known as the charge radius. 
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C is the convexity parameter given by 



C—-^r\y=i (28) 



The calculation of this IW function is non-perturbative in principle and is performed for different 
phenomenological wave functions of mesons [ 22 ] . This function depends upon the meson wave 
function and some kinematic factor, as given below : 

£(y)= Anr 2 \^(r)\ 2 cos(pr)dr (29) 
Jo 

where cos(pr) = 1 — + + with p 2 = 2p 2 (y — 1). Taking cos(pr) up to 0(r 4 ) we 

get, 



^ 2 2 4 4 

POO 'D ? 71 T* 

m = J 47rr 2 |*(r)| 2 (l-^ + ^)dr(30) 

= I" 4nr 2 \V(r)\ 2 dr - f°° ^r 2 \^{r)\ 2 [2p 2 ^ ~ ^- ]dr + [°° ^r 2 \^(r)\ 2 [Ap ^ V ~ ^ V ]rfr(31) 
jo jo 2 Jo 24 

= / Anr 2 \m(r)\ 2 dr - [Anp 2 / r 4 |^(r)| 2 dr](y - 1) + [-tt// 4 / r 6 |^(r)| 2 dr](i/ - 1) 2 (32) 
JO JO ' 3 Jo 

Equations (26) and (32) give us : 

roc 

p 2 = [Anp 2 r 4 \q{r)\ 2 dr] (33) 
Jo 

2 r°° 

C=[-Tip 4 r 6 \^(r)\ 2 dr] and (34) 

3 jo 



roc 

/ AKr 2 \m(r)\ 2 dr = 1 (35) 
Jo 



Equation (35) gives the normalization constants N and A/"' for \P°(r) and ^^(r) respectively 
, as : 



N (Jo 00 MiPir + Po])V 2 (36) 

N/ = [^/o°° r 2 (^A,( Pl r + Po ) - ^f(f + + 5p (37) 

3 Calculation and result: 

With linear confinement term of potential as parent, the wave functions contain Airy's function 
Ai[p±r + po] , which is an infinite series in itself [24]. 

A[ Pl r + Po ] = a[l + + __ + __ + ...]_ 

b[( Pl r + p ) + (Pir + Po)4 + iPir + Po)7 + (Pir + Po)1 ° + ...] (38) 

w 12 504 45360 J V ; 
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with a = 0.3550281, b = 0.2588194. 

Here we have studied the sensitivity of the order of polynomial approximation of the Airy's 
infinite series taking polynomial orders r 3 ; r 4 ? r 6 ; r 7 ? r 9 ; r 10 ( as polynomial orders r^ 2+3 ^ with 
1=0,1,2 etc are absent in the Airy's function series). Further, it is found that the infinite limit 
of integration in calculating and its derivatives makes the result divergent. Here, we take 
some reasonable cut-off limit tq of its integration. This will not sacrifice the nature and value of 
Airy's function and its derivatives, because, Airy's function falls very sharply ( exponentially) 
and almost dies out with increasing r- value. (Table 1 and Fig. 1,2 ). 



Table 1: Values of Airy's function for some small positive r. 



r 


Ai 


r 




r 


Ai [r] 


0.1 


0.329 


2.5 


0.016 


0.5 


0.232 


3.0 


0.007 


1.0 


0.135 


3.5 


0.003 


1.5 


0.072 


4.0 


0.0009 


2.0 


0.035 


4.5 


0.0003 



We find that Airy's function value becomes negligibly small for r > 5. Also, the graph of 
normalization constants (N and N' )versus the cut-off to upper limit ( r ) shows that N and 
N' value decreases with increase in r$ ( Fig. 3,4 ). 

Also, the graph of p 2 vs r ( Fig 5 ) and C vs r ( Fig 6 ) confirms that, beyond r = 9, p 2 
and C values rise steeply as compared to the result of Table 2. 



Table 2: Results of slope and curvature of £(y) in different models and collaborations. 



Model / collaboration 


Value of slope 


Value of curvature 


Ref [6] 


0.7936 


0.0008 


Le Youanc et al [25 ] 


> 0.75 


> 0.47 


Skryme Model [26 ] 


1.3 


0.85 


Neubert [27 ] 


0.82 0.09 




UK QCD Collab. [28] 


0.83 




CLEO [29,30 ] 


1.67 




BELLE [31 ] 


1.35 




HFAG [32 ] 


1.17 ±0.05 




Huang [33] 


1.35 ±0.12 





Upon this consideration, we have explored the £ (y) and its derivatives for different orders of 
polynomial approximation of Airy's function both for unperturbed wave function (Table 3) and 
total wave function (with relativistic effect) (Table 4), taking different cut-off values ranging 
from r = 5 to r = 9 , for D meson taking a s = 0.22 [23]. 

Regarding sensitivity of the order of polynomial in infinite Airy's function, the result for 
p 2 and C do not differ much upon variation of order of polynomial in Airy's function. For a 
given Airy order, with increase in cut-off value, p 2 and C values increases steadily, whereas for 
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Table 3: Result with unperturbed wave-function with r = 5,7,9 GeV 



0(r)in 




r = 5 






r = 7 






r = 9 




Ai 


N 


P 2 


C 


N 


P 2 


C 


N 


P 2 


C 




0.9307 


0.5053 


0.0872 


0.8864 


0.7595 


0.2213 


0.8656 


0.9857 


0.4581 


^4 


1.0658 


0.6712 


0.1123 


1.0121 


0.8625 


0.2195 


1.0118 


0.8658 


0.2229 


r 6 


1.0374 


0.6414 


0.1066 


0.9865 


0.8343 


0.2161 


0.9835 


0.8634 


0.2468 




1.0201 


0.6221 


0.1031 


0.9720 


0.8078 


0.2071 


0.9711 


0.8160 


0.2153 




1.0235 


0.6258 


0.1037 


0.9749 


0.8128 


0.2086 


0.9737 


0.8236 


0.2196 


r w 


1.0248 


0.6274 


0.1040 


0.9760 


0.8144 


0.2090 


0.9750 


0.8242 


0.2189 



Table 4: Result with total wave-function ( with relativistic effect) taking ro = 5,7,9 GeV 



0(r)in 




r = 5 






r = 7 






r = 9 




A 


N 


P 2 


C 


N 


2 

P 


C 


N 


2 

P 


C 




1.5927 


0.5149 


0.0661 


1.5658 


0.6125 


0.1269 


1.4985 


0.9681 


0.4365 


^4 


1.8653 


0.6300 


0.0869 


1.8104 


0.7942 


0.1932 


1.4375 


2.4989 


1.8879 


r 6 


1.8703 


0.6432 


0.0878 


1.8162 


0.8040 


0.1921 


1.5311 


2.0925 


1.4440 


r 7 


1.8471 


0.6313 


0.0858 


1.7946 


0.7899 


0.1884 


1.4769 


2.2604 


1.6355 




1.8518 


0.6340 


0.0862 


1.7989 


0.7931 


0.1892 


1.4871 


2.2294 


1.5990 


r w 


1.8531 


0.6348 


0.0863 


1.8000 


0.7940 


0.1895 


1.4847 


2.2462 


1.6169 



a given cut-off value, p 2 and C values do not differ much upon variation of order of polynomial 
in Airy's function from r 4 to r 10 . However, the results show closer resemblance to recent result 
of p 2 = 1.17 [32] for our Airy-order r 3 up to cut-off value r = 10. For such specific range and 
order, our result shows improvement upon the result of ref [6]. At cut-off value higher than 
r = 9, the results jumps to higher values than our expectations. In all calculations, we have 
taken care that the boundary condition of IWF ( £(1) = 1 ) is maintained throughout. The 
variation of £(y) with y for different cut-off value r° ,with Airy order r 10 , is shown in figure 
7 and the variation of £(y) with y for different Airy order at cut-off value r = 5 is shown in 
Fig. 8. In 8, graphs of vs y overlaps for Airy orders r 4 to r 10 , whereas the graph for Airy 
order r 3 shows a slight deviation. 



4 Conclusion and remarks: 

We have found that cutting off the upper limit of integrations in £(y) and its derivatives to 
some reasonable point does not upset the result, rather it almost conforms to the experimental 
expectations. Also, for each value of cut-off r , we have considered the asymptotic form of the 
Airy's function taking limits of integration from r to oo . 

ex P (~| r3/2 ) 
Ai[r\ asympt ~ 2 v /7rr 1 / 4 

With this asymptotic form we have also calculated the derivatives of £(y). Such analysis shows 
that very small values of p 2 and C result [Table 5], taking this asymptotic form of Airy's func- 
tion. This, thus, also confirms the justification of taking some reasonable cut-off to upper limit 
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of integration in calculating £ (y) and its derivatives. 



Table 5: Values of p 2 and C with asymptotic form of Airy's function. 



r value 


p 2 (asymptotic ) 


C (asymptotic) 


5 


4.6 x 10~ 9 


1.6 x 10~ 9 


6 


5.027 x 10~ n 


2.464 x 10~ n 


7 


3.56 x 10" 13 


2.345 x 10" 13 


8 


1.695 x 10" 15 


7.028 x 10" 15 


9 


5.248 x 10" 18 


6.597 x 10" 15 


10 


2.92 x 10" 19 


2.78 x 10" 19 



Plots of versus y for different orders of Airy's function are also in agreement with our 
expectations [8]. The graph of with y for different Airy's function order invariably starts 
at (1,1) and almost follows the same pattern and show very small deviation with change in Airy 
order. It confirms the fact that boundary condition for zero recoil (£(1) = 1) is maintained all 
through, with different polynomial orders of Airy's function and for different cut-off values. 
Let us also comment on the result of Ref [6]. The result of ref [6], which is for Airy order r 3 , 
matches with our cut-off value tq = 7.95 GeV^ 1 for the same Airy order in our calculation. 
However, their wave function is not found to satisfy the zero recoil condition of IW function, 
f (y) = 1 for y = 1. 

To conclude, we also study the compatibility our potential model with the recent results of 
Heavy Flavour Averaging Group [32]. Taking the result of p 2 , we fix the value of cut-off for 
different orders of polynomial in Airy function (Table 6). It indicates that, the range of cut-off 
value r = 7 GeV^ 1 to r = 9 GeV^ 1 matches the expectations of ref [32]. 



Table 6: Value of cut-off r for different Airy order matching expectation of p 2 . 



p 2 = l.nGeV- 1 (ref [32]) 


Airy Order 


r value 




8.896 


r 4 


7.915 


r 6 


7.975 


r 7 


7.942 


r 9 


7.939 


r io 


7.932 



Lastly, we would also like to comment on the limitation of the present formalism. The per- 
turbed wave function \l/ 1 (r) has been calculated using Airy order r 3 , as in ref [6]. Although, the 
total wave function (eqn 22) contains infinite Airy series in terms of unperturbed wave function, 
this above mentioned limitation may have some effect on the values of the coefficients of the 
higher polynomial orders of Airy function in total wave function. Improvement of the formalism 
considering higher polynomial orders of Airy's function in ^(r) is under consideration. 
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Figure 1: Variation of Airy's function with r 
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Figure 2: Variation of Airy's function with positive r 
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N vs rO 
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Figure 3: Variation of N with r . 



N 1 vs rO 
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Figure 4: Variation of V with r . 

Slope v/s cut-off with Airy-r^lO 



. 6 
^0.4 
0.2 



2 4 6 8 

rO 

Figure 5: Variation of p 2 vs r with Airy order r 10 . 
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Curvature v/s cut-off with Airy ~r A 10 
0.1 




Figure 6: Variation of C vs r with Airy order riO. 




Figure 7: Variation of £(y) with y taking different r with Airy order r 10 . 
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Figure 8: Variation of with y taking different Airy order with r = 5. 
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